Pressure dependence of critical temperature in MgB2 and two bands Eliashberg theory 
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The variation of the superconducting critical temperature Tc as a function of the pressure p in the 
magnesium diboride MgB2 has been studied in the framework of two-bands Ehashberg theory and 
traditional phonon coupling mechanism. I have solved the two-bands Eliashberg equations using 
first-principle calculations or simple assumptions for the variation, with the pressure, of the relevant 
physical quantities. I have found that the experimental Tc versus p curve can be fitted very well and 
information can be obtained on the dependence of the electron-phonon interaction matrix < I^ > 
by pressure. The pressure dependence of the superconductive gaps and is also predicted. 

PACS numbers: 74.62.Fj; 74.62.-c; 74.20.Fg; 74.70. Ad, 74.25.Kc 
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In the last few years, there is an noticeable increase 
of the study of superconductivity in many elements un- 
der pressure 0, such as sulphur {Tc = 17 K), oxygen 
(Tc = 0.5 K), carbon in nanotubc (Tc — 15 K) and di- 
amond forms (Tc = 4) K, a non-magnetic state of iron 
(Tc = 1 K), and the light elements lithium (Tc = 20 
K) and boron (Tc = 11 K). The application of ex- 
ternal pressure to superconductors can drive the com- 
pounds towards or away from lattice instabilities by vary- 
ing the principal parameters determining the supercon- 
ducting properties (the electronic density of states at the 
Fermi energy, the characteristic phonon frequency, and 
the electron-phonon coupling constant), and it can be 
used to tune the Tc and the superconducting properties. 
Almost all of the superconducting metallic materials, un- 
like the previous simple elements, show a decrease of Tc 
with pressure. This negative pressure coefficient was at- 
tributed to the volume dependence of the electronic den- 
sity of states at the Fermi energy and of the effective in- 
teraction between the electrons mediated by the electron- 
phonon coupling. Measurements of the influence of pres- 
sure on the transition temperature and critical field yield 
information on the interaction causing the superconduc- 
tivity. Indeed, the pressure would seem to be a variables 
whose effects might be capable of immediate theoretical 
interpretation. The binary alloy MgB2, superconductor 
0, at ambient pressure, at T = 40 K has, under pressure, 
a behaviour similar to metallic materials. The magne- 
sium diboride has stimulated intense investigation, both 
from the theoretical and the experimental point of view. 
Now the electronic structure of MgB2 is well understood 
and the Fermi surface consists of two three-dimensional 
sheets, from the tt bonding and antibonding bands, and 
two nearly cylindrical sheets from the two-dimensional 
a bands There is a large difference in the electron- 
phonon coupling on different Fermi surface sheets and 
this fact leads to a multiband description of supercon- 
ductivity. Theory indicates that the strongest coupling is 
realized for the near-zone center in-plane optical phonon 



{E2g symmetry) related to vibration of the B atoms Q]. 
The superconductivity in MgB2 has been deeply studied 
in the past three years and so also the effect of pressure on 
the superconductive properties. The effect of pressure on 
the superconducting properties of MgB2 has been stud- 
ied by several groups. All groups observed a decrease of 
Tc with increasing pressure 0,13 and I want show that 
this decrease can be very well explained in the framework 
of the two bands Eliashberg theory. In the following I will 
refer to the paper of A.F.Goncharov because in there 
are present both measurement of the variation of critical 
temperature and of phonon mode by means of Raman 
measurement, with the pressure and so I mainly refer to 
these experimental data. In fact only in this work there 
are all input parameters necessary to my model. 

Let us start from the generalization of the Eliashberg 
theory 0, |1] for systems with two bands |^, that has 
already been used with success to study the MgB2 and 
related systems (ItX ll. 12. 13. 14. 15j. To obtain the gaps 
and the critical temperature within the s-wave, two-band 
Eliashberg model one has to solve four coupled integral 
equations for the gaps Ai(ia;„) and the renormalization 
functions Zi{iujn): 



(1) 



Zi(ia;„)Ai(iw„ 



7rT^[Ay(ia;„ - iujrn) - fJ'iji'^c)] ■ 
) + J2^^Nii^co„) (2) 



where i,j are band indices, 9 is the Heaviside function, 
ujc is a cutoff energy, F*^ is the non-magnetic impurity 
scattering rate in the Born approximation and: 
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FIG. 1: The spectral functions of the two-band model for the 
MgB2: aa (solid line), nn (dashed line), an (dotted line) and 
na (dashed dotted line), taken from ref. 11. 



FIG. 2; Experimental energy of phonon E2g of MgB2 as a 
function of the pressure, see ref. 5. In the insert the exper- 
imental Raman linewidth as a function of the pressure, see 
always ref. 5. 
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where Wn = nT{2n — 1) and n, m = 0, ±1, ±2.... 

The solution of Eqs. 1,2 requires as input: i) the four 
(but only three independent electron-phonon spec- 
tral functions afj(w)F(w); ii) the four (but only three 
independent 0) elements of the Coulomb pseudopoten- 
tial matrix pL*{uJc)- 

Let's start with the four spectral functions af^ {lo)F{lo), 
that were calculated in ref. 11 (sec Fig. 1). 

For simplicity, I will assume that the shape of the 
afjF{u!,p) functions does not change with the pressure, 
and I will only rescale them with the electron-phonon 
coupling constants Ay : 



lJp)=N\{p)<f>IMn^ 'ip) 
xjp)=cm-^{p)a^y)" - 
X (p)=this model "' 
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FIG. 3: Calculated electron-phonon coupling constant as a 
function of pressure in the four different cases: Ao-ct(p) = 
n/^l^M (dotted line), \aa{p) = N^{Ev,p) < l'^ > 



/2Am%^^ip) (dashed line), A..(p) = C/(iV5(i5F,p)n|,^ (p)) 
(short dashed line) and \crcr{p) = [l{p)/^%2 (P) 
■kX^^{0)N^{Ef,0) - 'yiO)/nl^(0)]/7vN^{Ep,p) (solid" line). 
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Let me remind the definition of electron-phonon coupling 
constant 0,0, [iSj: 



(7) 



a band, from eq. 1 will have for A^ct (which gives the 
most important contribution to superconductivity in our 
system): 



Acrcr(p) 
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where "fi{q) is the phonon linewidth which is the width 
in energy of a phonon of momentum q, mode index i and 
energy fli{q) and N-^{Ep) is the normal density of states 
at the Fermi level. The frequency f2i can be identified 
with the frequency of the B-B bond-stretching phonon 
mode (the £^2g mode), that has been recently measured 
as a function of pressure 0, and is reported in Fig. 2. In 
the insert of Fig. 2 the experimental Raman linewidth 
that here is used as phonon linewidth, is shown. Since 
this mode couples strongly with the holes on top of the 
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where 



C{0)=7rX,,iO)N^{EF,0) 
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When the lattice parameters of MgB2 are modified by 
chemical substitutions the normal density of states, at 
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the Fermi level, in the 7r-band changes relatively little 
[l9i | and so I assume that, in the first approximation. 



and 



N^{EF,p)^N^iEF,0)+p 



dN^jEp.p)^ 
dp 



\p=o 



(11) 



(12) 



I use the values calculated in ref. 10: N^{Ef,0) = 
0.30061 {eVunitcell)-^ and N^{Ef,0) = 0.40359 



ajV5(gF,p)i 



=n is the 



lp=0 



{eV unitcell)~^ for the MgB2. So 
only true free parameter of the model. In this way, I as- 
sume that the change in the frequency of the i?2g phonon 
affects the coupling constant, while I neglect its influence 
on the shape of the electron-phonon spectral function 
[T^ . For the other coupling constants, I will instead as- 
sume for simplicity 



V(«, j) ^ (o-,cr) 



A^(i^P,p) 



(13) 



with [13, [T3 A<,<,(0) = 1.017, A^^(0)=0.448, 
A^^(0)=0.213 and A^^(0)=0.155. At the end, in this ap- 
proximate model of electron-phonon coupling constants 
only Aero- and A^o- change with the pressure. This fact is 
in agreement with the results of ref. 20 where the au- 
thors find that Atttt is almost constant. Fig. 3 shows 
the calculated electron-phonon coupling constant X^a as 
a function of the pressure. 

As far as the Coulomb pseudopotential is concerned, 
let us start from its expression in pure MgB2 (10t.»lL . i2]| : 



2.23 



M(c^c)A^r(^F,p) 



1 



1 2.48 



(14) 



where ^{coc) is a free parameter and N^'^{E-p,p) is the 
total normal density of states at the Fermi level. The 
numbers 2.23 and 2.48 in the Coulomb matrix have been 
calculated for the MgB2 in ambient pressure but, as a 
first approximation, I will suppose them not to depend 
on the pressure. In this way, the elements of the Coulomb 
pseudopotential matrix, //jj , depend on the pressure only 
through the densities of states at the Fermi level. Now 
I can fix the cut-off energy (e.g., uJc ~ 700 meV) so as 
to reduce the number of adjustable parameters to two: 
the prefactor in the Coulomb pseudopotential, ii{lOc) and 



dN^{EF,p) 
dp 



\p=o- For having Tc — 40.2 K I fix /^(wc) equal 
to 0.00315 so there is only more a free parameter for fit- 
ting the experimental critical temperature as a function 
of pressure. Before of examining the fit of experimental 
data with the my Ao-o- model I can check the other pos- 
sible choices for the electron-phonon coupling constants. 
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FIG. 4: The measured critical temperature Tc as func- 
tion of the pressure (open circles) and theoretical fits ob- 
tained by different assumptions on the electron-phonon 
coupling constant Xa-a-ip) 



(p) (dotted line). 
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j{0)/Ql.^m/nNUEF,p) (solid line). 



All cases are shown in Fig. 3. The first and simpler pos- 
sibility is similar to the case of several transition metals 



X{p) = 



(15) 



where r/ is a constant [1^ |23 and so 



Xaa{p) = 



"k(0) 



(16) 



The result is in very poor agreement with experimental 
data (see dotted line in Fig. 4). 

The second possibility is of that the effect of the pres- 
sure is similar to chemical substitutions as Al and C and 
so I assume that [Tj 



Kp) 



n^{Ef,p)<P> 
Mnl (p) 



(17) 



where Af is the ion mass and < P > docs not depend 
from the pressure. Consequently 



Xacr{p) 



N-{ET.,p)nljQ) 
N-{ET.,0)nl (p) 



(18) 



ajv5gF,p) |^^^ ^ -0.003 {eVGPa)-\ Also in this case 
the result is in very poor agreement with experimental 
data (see dashed line in Fig. 4). 
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FIG. 5: The calculated values of the gap, at T = 4.2 K, A<j 
(solid line) and (dashed line) as a function of the pressure. 



The last possibility is suggested by recent band- 
structure calculations that show MgB2 is a traditional 
sp metal superconductor The pressure dependence 
of / has long been an interesting issue in the research of 
pressure effects in simple sp metals |2^. Zimans calcu- 
lation of the electron-phonon interaction led to < >cx 

7Vn(Sf)-~ 
So I find 



at least in the limit of long wavelengths [22 



jVg(gF,0)»|,^(0) 
N-{E^,p)al ip) 



(19) 



and 



3Afg(i;F,p)i 



-g^ |p=o = -0.0007 [eVGPa)-^. As in the 

previous cases the result is in very poor agreement with 
experimental data (see short dashed line in Fig. 4). 
Now I can see that my simple model is the alone that 
explains the experimental critical temperatures because 
other possible models for electron-phonon coupling con- 
stant A(j(j are incompatible with experimental data. I 
obtain the best fit of experimental data (solid line in Fig. 



4) with 



3jVg(£;F,p)| 



0.00584 {eVGPa)-^. The fact 



that N^{E-p) increases with t he p ressure is in agreement 



nepr 

with theoretical calculations [22|. The only free param- 



eter of this model is 



p=o and so when I have 



fixed the optimal value from the Tc fit I can calculate, in 
principle, all other physical quantities. In Fig. 5 the the- 
oretical dependence of the a and tt gaps from the pressure 
is shown. 

Now from the following equality 



N^iEF,p) <P{p) > 



1 



r7E2g(p) 



nN-{EF,p)'nE, 



+ 



nil)- 



(20) 

it is possible determine the dependence of N-^{E-p) < 
P > and of < > from the pressure as it is shown 
in Fig 6. It can see that, after ~ 25 GPa the Hopfield 
parameter |23| N^{Ef) < > is almost constant. 

At the end it is possible to use this model for explain- 
ing the experimental upper critical field '20i| in function 
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FIG. 6: Upper panel: the calculated value of N^{Ef) < > 
of the (T-band as a function of pressure; lower panel: the calcu- 
lated value of < > of the a-band as a function of pressure. 
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FIG. 7: The measured upper critical field (open circles), at 
T=10 K, from ref. 20, and the theoretical (solid line). 



of pressure without free parameters. For the sake of com- 
pleteness, I give here the linearized gap equations under 
magnetic field, for a superconductor in the clean limit 
(negligible impurity scattering), as can be found in ref. 
20. In the following, vpj is the Fermi velocity of band j, 
and Hc2 is the upper critical field: 

uJnZi{iuJn) = t^ri + ttT ^ Ay (iw„ - iujm)sign{ujrn) 

(21) 
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Zi{iu!n)Ai{iuJri) = ttT ^ [Ay (iw„ - it^™) - /X,*,- (wc)] • 
■6{\uJc\ - UJm)Xj{i^m)Zj{iuJm)Aj{ii^m) 

(22) 

r + co 

Xj{i^m) = (2/\/^) / d(jexp(-q^)- 

-'0 




\ujmZj{iu;m)\ +ifJ.BHc2sign{u!m) 

(23) 

with f3j — 'KHc2Vpj/{2^o). In these equations the 
bare Fermi velocities are the input parameters vpj — 
Vpj ■ + ^ji) ^iid are functions of p. For don't 

having free parameters I assume that, as in the 
free electron gas, vpj cc N^{Ef) and so VFj{p) = 
VFj{0)N^iEF,p)/N^{EF,0). For obtaining exactly the 
upper critical field of MgB2, in ambient pressure, I find 



v*p„{Q) = 3.6 ■ 10^ m/s and v*p^{Q) = 5.35 ■ 10^ m/s in 
very good agreement with the calculus of Brinkman et al 
. In Fig. 7 is shown the fit (solid line) of experimental 
values, from ref. 26, of Hc2 at T = 10 K. The fit isn't 
so good because, may be, the approximation of the free 
electron gas is too strong. 

Finally I conclude by summarizing the main points of 
this paper. I have fitted the experimental critical temper- 
atures as a function of pressure in the framework of two 
bands Eliashberg theory with only a free parameter. The 
result is very good. After I have calculated other physi- 
cal quantities can will be compared with future measure- 
ment (for example superconductive gaps from tunneling 
curves) and I can affirm that the MgB2 under pressure 
is, as the same materials in ambient pressure, a moder- 
ate coupling two-band phononic systems well described 
by two-bands Eliashberg theory. 

Many thanks are due to A. Calzolari for useful discus- 
sions. 
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